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The spin-orbit interaction affects the electronic structure of solids in various 
ways. Topological insulators are one example where the spin-orbit interaction 
leads the bulk bands to have a non-trivial topology, observable as gapless surface 
or edge states. Another example is the Rashba effect, which lifts the electron- 
spin degeneracy as a consequence of spin-orbit interaction under broken inver- 
sion symmetry. It is of particular importance to know how these two effects, i.e. 
the non-trivial topology of electronic states and Rashba spin splitting, interplay 
with each other. Here we show, through sophisticated first-principles calcula- 
tions, that BiTel, a giant bulk Rashba semiconductor, turns into a topological 
insulator under a reasonable pressure. This material is shown to exhibit several 
unique features such as, a highly pressure-tunable giant Rashba spin splitting, 
an unusual pressure-induced quantum phase transition, and more importantly 
the formation of strikingly different Dirac surface states at opposite sides of the 
material. 

Introduction 

Theoretical work of Kane and Mele [lj] has marked a milestone in our understanding of the 
insulating phase of the matter. In their seminal work, they classified the bulk insulators 
according to the topological order of their electronic states using quantum metrics known 
as topological indices. When nontrivially ordered, it was predicted that such topological 
insulators (TPs) should exhibit gapless states at their boundaries, e.g. at the surface or 
edges. The succeeding works have confirmed this prediction by theoretically proposing 
and/or experimentally discovering a number of two and three dimensional TI's, including 
HgTe quantum wells 

as 

and uniaxially strained bulk HgTe , a number of layered 



and the ternary compounds 



urt 



ner 



V2VI3 binary compounds e.g. Bi 2 Se3, Bi 2 Te3, Sb2Te3 {(J 

LaBiTe3 TlBiTe2, and TlBiSe2 9, 10 1. The recent first-principles studies have 
proposed that several half-Heusler compounds such as LaPtBi, LuPtSb, and YPdBi [H 
can turn into a TI under a uniaxial strain, similar to what was predicted and observed 
for HgTe % l5|- Common to all these materials is the presence of heavy elements with a 
reasonably large atomic spin-orbit interaction (SOI). Such a similarity arises from the fact 
that in all TI systems the bulk band gap is the result of band inversion by SOI [4]. In other 
words, without such a SOI-induced band inversion, no topological phase can be realized in 
a system. 



2 



In a more general context, SOI tends to lift the degeneracy of energy bands. In systems 
with inversion (J) symmetry, such a splitting leaves the energy states spin-degenerate, as 
long as the time-reversal (T) symmetry is hold. T-symmetry connects the states an d 
V'-k,^) while the /-symmetry enforces the degeneracy between ^k.td) an d V'-k,tGQ- Breaking / 
symmetry lifts the latter constraint and, hence, lets the energy bands be spin-split at generic 
fc-points. Rashba spin splitting (RSS) 15[ well exemplifies this situation, as described by 
Hr = 2— + z/[e ■ (s x p)] where e is the direction of the potential gradient, breaking I- 
symmetry, and s and p are the spin and momentum operators, respectively. This interaction 



leads to several unique phenomena, such as the spin Hall effect 16J, the spin Galvanic 
effect 17[, and the magneto-electric effect 18[. Furthermore, once superconductivity occurs 



in a Rashba system, unusual features such as the mixing of the singlet and triplet pair states, 



a large upper critical field beyond the Pauli limit [19|, [20| , and topological superconductivity 
with Majorana edge channels 2l|] can appear. As SOI is responsible for both RSS and 
non-trivial topology of electronic states in semiconductors, it is therefore of fundamental 
interest to know how they interplay with each other once they coexist in a system and more 
importantly what the impact of such an interplay is on the electronic structure of bulk and 
boundary states. 



The purpose of this work is to study such an interp 



ay in a layered polar compound, BiTel. 



Backed by our earlier band structure calculations [22|, |23[ , the angle resolved photoemission 
spectroscopy (ARPES) measurements 22, [24] have revealed that the bulk conduction and 
valence states in BiTel are subject to a giant RSS of the order of several hundred meV, 
lying among the highest reported so far. In this work we further show using the first- 
principles calculations that by applying a reasonable hydrostatic pressure, P, the material 
turns into a TI with many interesting features, such as a nearly double enhancement in RSS 
accompanied with an unconventional metallic behavior at quantum phase transition. It is 
also demonstrated that unlike the /-symmetric TPs, the gapless Dirac states in TI phase 
of BiTel have completely different shapes at different sides of the material. Consequently, 
here the surface states on both sides can interestingly have the same spin helicity in a rather 
broad range of energies inside the bulk band gap, in sharp contrast with the centrosymmetric 
TPs. 

Results 



Basic properties. BiTel belongs to the trigonal space group of P3ml. As shown in Figs. 
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l-(a) and l-(b), the crystal structure of BiTel has a non-centrosymmetric layered structure 
along its crystallographic c-axis with three atoms in one unit cell. Within each unit, a Bi 
atom is sandwiched between one Te and one I, forming a triple layer. Due to the strong 
covalency and ionicity of Bi-Te and Bi-I bonds, respectively, the bulk crystal intrinsically 
possesses a polar axis along the ^-direction. Despite the strong chemical bonding within each 
triple layer, the adjacent triple layers are weakly coupled via van-der Waals interaction. From 



our previous calculations 



231 ] . we know that around the Fermi level Ep, all the bands are 



essentially p-type, with conduction bands dominated by Bi-Qp and valence bands composed 
of Te-5p and I-5p states. Moreover, due to the negative crystal field splitting (CFS) of the 
valence bands and positive CFS of conduction bands near Ep, in the absence of SOI the 
top valence bands (TVB's) and bottom conduction bands (BCB's) both become p z type 
(to avoid any confusion, hereafter they are referred to as p A and pf , respectively). These 
features make the electronic structure of BiTel very much similar to that of the well known 
TI systems Bi2Te3 and Bi2Se3 However, here due to the absence of /-symmetry, it's not 
possible to assign a distinct even or odd parity to each band. It is also worth noting that 
in the case of BiTel, the minimum energy gap Eg is not at the Brillouin zone (BZ) center 
T but areound the hexagonal face center of BZ, A point, where k x = k y = and k z = tt/c, 
as shown in Fig. l-(c). As will be shown below, because of the latter difference, BiTel in 
TI phase shows a rather different topological order from what has been found for Bi 2 Te3 
and Bi2Se3. Introducing SOI, both spin and orbital mixings are allowed. Consequently, 
pf (pf) transforms to \p A , ±|) ( |p B ,±|)), thereby getting energetically repelled upward 
(downward). This accordingly closes the band gap from 1.2 eV down to 0.286 eV and 
induces a giant bulk RSS among these two sets of bands around point A (see discussion in 



Ref. 



23] ). Despite such a huge reduction in Eg, BiTel remains a trivial insulator as Eg still 



originates from atomic orbital hybridization between Bi and its neighboring Te and I atoms. 
The respective band diagram is shown in Fig. l-(d). 

Topological phase transition in BiTel under pressure. Our strategy to turn BiTel 
into a TI is to modify its chemical bonds by applying an external hydrostatic pressure. 
As schematically shown in Figs. l-(d), l-(e) and l-(f), through this modification we can 
effectively control both CFS and SOI such that at a critical pressure P c , \p A , ±|) and \p B , ±|) 
become degenerate, whereas at higher P's their energy ordering is reversed thereby forming 
an inverted band gap. Controlling CFS by P is rather easy to understand, because any 
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change in Bi-Te and Bi-I bonds leads to a change in the energy splitting of pf ,B and p^'y 
states. For example, our non-relativistic band structure calculations reveal that at P c , the 
CFS of TVB's is so enhanced that, Eq is reduced by 200 meV (see the Supplementary Fig. 
SI). Any band gap narrowing associated with CFS modification can be further enhanced 
through a subsequent enhancement of RSS of \p A , ±|) and \p B ,±~) states. As described 



in detail in Ref. 



23] , this is due to the fact that these two states are symmetrically of the 



same character and hence can very effectively couple with each other through a Rashba-type 
Hamiltonian if they are energetically close to each other. In other words, the closer they are 
to each other, the larger RSS would be achieved. 

To elucidate this mechanism, we show in Figs. 2- (a), 2-(b) and 2-(c) the respective 
electronic band dispersions of TVB's and BCB's along the high symmetry direction H—A—L 
of BiTel as hydrostatically compressed by V/Vq = 1, V/Vq = 0.89 and V/Vq = 0.86, where 
V corresponds to the lattice volume at the ambient pressure P am bient- As shown, at P am uent 
(V/Vq = 1) a comparable giant RSS can be seen for both sets of bands with an Eg of 286 
meV. The corresponding Rashba energy Er, defined as the energy difference between the 
conduction band minimum (CBM) and the conduction band crossing point, is nearly 110 
meV, in perfect agreement with that observed by ARPES measurements j^j. Compressing 
the volume down to V/Vq = 0.89, Er monotonically increases until it reaches 200 meV, 
astonishingly about two times larger than that at Pambient, as shown in Fig. 2-(b) (for 
a detailed comparison see the Supplementary Fig. S2, also). At this point, the system 
reaches its quantum critical point, represented by a full band gap closing along the A — H 
directions. It is worth noting that such a quantum phase transition in BiTel differs from 
that in usual centrosymmetic TI's, as there it is mediated through a band gap closing at 
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2611 . whereas in BiTel as will be shown 



a single high symmetry /c-point, e.g. T point 
shortly due to the spin splitting the band gap is closed at six Appoints, all along A — H 
directions. Thus, they are describable by six 2 x 2 Dirac Hamiltonians. Here, it is important 
to stress that because of T-symmetry these six fc-points are two-by-two paired with each 
other, therefore their corresponding states form three pairs of 2 x 2 Dirac fermions. Since 
the band touching is mediated through odd number of such pairs, the transition to TI phase 
is in principle allowed. Further increasing P, Er starts decreasing and an inverted band 
gap emerges between TVB's and BCB's. Figures 2-(d), 2-(e) and 2-(f) demonstrates the 
orbital characteristics of BCB's and TVB's at V/Vq = 1, V/Vq = 0.89 and V/Vq = 0.86, 
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respectively. A thorough analysis of these bands reveals a clear change in their atomic 
orbital characteristics as a function of V/Vq. While at V/Vq = 1 the BCB's and TVB's are 
dominated by Bi-6p and Te-5p, respectively, at V/Vq = 0.89 they appear to share rather the 
same orbital characters near the band touching point. On the other hand at V/Vq = 0.86, 
one can clearly see that Bi-6p orbitals now contribute more significantly to TVB's whereas 
the BCB's become strongly dominated by Te-5p and I-5p states, a clear indication that a 
band inversion has happened and, hence, the system is now in the TI phase. It is worth 
adding that in other parts of BZ including the plane encompassing T, M and K the band 
gap remains finite for the whole range of P's. 

To make sure that BiTel for P > P c becomes a TI, we have calculated the so-called Z 2 
topological invariant of the whole system when hydrostatically compressed by V/Vq = 0.86. 
By definition, for three-dimensional systems, Z 2 topological invariant is determined by four 
indices, among which one is called the strong topological index v and the three others 
are weak topological indices u±, v-i and 1/3 [4j. The complete list is then described as 
Z2 = Vo] (^1.^2^3) • For a normal insulator all the four indices are zero whereas a strong 
TI is characterized by vq = 1. If vq = but any of ^1-3 is 1, the corresponding system is 
called a weak TI. For the numerical computation of Z 2 invariants, several approaches have 



been proposed [27H30|. Here we follow the approach introduced in Ref [271 (see the Supple- 



mentary Methods and Supplementary Fig. S3 for a full description of our computational 
method for the calculation of the Z2 topological invariants). Our calculations show that at 
V/Vq = 0.86, BiTel belongs to Z2=l;(001). This clearly indicates that it is a strong TI. 
Interestingly, the corresponding ^3 turns out to be 1, which can be attributed to the band 
touching in the hexagonal face of BZ with k z = tt/c, encompassing A and H points. The 
interference between the physics of strong and weak topological insulators could be there- 
fore an interesting problem here. For example, V3 = 1 means that the one-dimensional (ID) 
conduction channel can appear by dislocation with a Burger's vector along the z-axis 31]. 
Thus, the present system may offer an ideal medium to study the interaction between the 
surface Dirac fermion as guaranteed by vq = 1 and such a ID channel. 

To further confirm that BiTel becomes a TI for P > P c , we next show in Fig. 3 the 
surface band structures corresponding to I- and Te-terminated sides of BiTel. While the 
surface states are fully gaped at P a mMent, a gapless state appears for the both sides at P c . 
At higher pressures, as shown for the case of V/Vq = 0.86, gapless surface states appear 
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within the bulk band gap. A clear indication that BiTel has now become a strong TI 
and, thereby further confirming the results of our Z2 invariants calculations. Interestingly, 
the shape of the Dirac surface states at the Te-ended side of BiTel completely differs from 
that at the I-ended one. While for the former the Dirac point is deeply buried inside the 
energy valley formed by the Rashba-split TVB's (resembling the Dirac surface states in 
Bi2Te3), on the other side the Dirac point is energetically well above the bulk CBM. By 
increasing P, due to the reduction of RSS and the subsequent increase in Eg, the latter 
Dirac point becomes energetically closer to CBM. We expect this to be a characteristic 
feature for all non-centrosymmetric TI candidates making them distinguishable from I- 
symmetric TI systems. As a consequence, in a non-centrosymmetric TI, the electrons are 
fractionalized into nonequivalent halves on top and bottom surfaces. This would lead to 
some novel features and new spintronics functionalities. As shown in Fig. 3-(g) and 3- 
(h), the surface Dirac fermions show the similar spin polarization patterns at the top and 
bottom surfaces, in sharp contrast to the centrosymmetric TI's. It is to be noted that for 
some thin slabs of Bi2Se3 and Bi2Te3, such an effect appears to happ en when the equivalence 
of the two surfaces is lifted, e.g. by applying an electric field 32[ or by differing chemical 



functionalization 



33(. Because of this effect, the Dirac point traverses across the bulk band 



gap on going from the top side to the bottom side, and always crosses Ep on the side surface 
as long as Ep is within the bulk gap. Therefore, if one applies a magnetic field or dopes 
magnetic impurities, an insulating stripe is expected to appear on the side surface (since Ep 
is now located inside the energy gap of the side surface Dirac fermion). Also the giant spin 
Galvanic effect is expected since the current-spin relation of the top side is the same as taht 
of the bottom side. When the hybridization between the Dirac fermions at top and bottom 
surfaces occurs, in the momentum space a gap appears along a one-dimensional (nearly) 
circular path, resulting in the divergence of the density of states. This accordingly leads 
to the enhanced electron correlation effect and, consequently, to the excitonic instability. 
These are just a few examples of the novel phenomena expected in the non-centrosymmetric 
TI's. 

Estimation of P c . Let us now briefly discuss on the possible value of P c required for 
topological phase transition. Performing two sets of volume optimization calculations using 
the local density approximation (LDA) and generalized gradient approximation (GGA) and 
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then fitting the respective free energies E(V) to the Murnaghan equation of state [34] 



E(V) = E + 



B V 
B' 



{Vo/V) 



B' -l 



+ 1 



B Vo 



B'-l 



(1) 



(see Fig. 4), we estimate the upper and lower limits of Bulk modulus Bq of BiTel to be 
~ 21.9 GPa and ~ 8.9 GPa. The respective LDA and GGA values of pressure derivative of 
bulk modulus B' Q are similarly found to be 7.7 and 7.8. As mentioned above, our calculations 
indicate that at P c , V is compressed by 11%. Using the relation 34] 

Bo 



P(V) 



B'o 



Vo} B '° 
V 



(2) 



P c is accordingly expected to be in the range of 1.7 GPa to 4.1 GPa. Given that this range 
of P's is rather easily affordable in laboratory, we thus hope this work would stimulate 
experimentalists in this field to explore such an intriguing topological insulating phase in 
BiTel. 

Discussion 

At this point it is worth explaining (i) as to why BiTel becomes gapless only at a certain 



P c but not for a range of P's as discussed before |35|, |36| and (ii) why such a gapless state 
occurs along a specific direction. Group theory is the key to answer these questions. Let 
us first explain the role of crystalline symmetry for the development of dispersion minima 
(maxima) at BCB (TVB) along particular directions. As already mentioned, BiTel belongs 
to C% v symmetry consisting of a three-fold rotation C3 along the ^-direction and three mirror 
operations M: y — > —y where y is along A — H directions. For spin 1/2 electrons, C3 and 
M can be represented as e -10 "- 71 "/ 3 and ia y , respectively, where cr x ,y,z are Pauli matrices for 
spin degrees of freedom. Additionally T operator can be defined as ia y K, where K is the 
complex conjugation. We can then construct a two band Hamiltonian H c (k) for the BCB's 



invariant under C3, M and T. Up to cubic terms of k it turns out to be: 



H c (k) 



K + k 2 y , (k-TT/c) 



2m 



+ 



_L,c 



2ml 



+ Vk,c( k x a y ~ K a ^i + ~ k 2 y )k y a z , 



(3) 



where m* ± c and m| c are the in-plane and out-of-plane effective masses of BCB's and = 
v c {\ + a c k 2 ). Note that in a similar manner, one can construct H v (k) for TVB's. The third 
term in H c (k) is obviously the Rashba term allowing a cylindrical in-plane spin splitting 
within k X;V plane. The fourth term, hereafter referred to as H w (k), acts as a warping term, 
trigonally distorting the energy bands. 
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Due to H w (k), the inner and outer branches of the Rashba-split conduction bands are 
distorted in a way that by approaching CBM they first merge together at fc-points along 
A — L directions and then form six energy pockets each of which centered along one of six 
A — H directions, as shown in Fig. 5- (a). Exactly at CBM these six pockets reduce to six 
points. In the same manner, H w (k) distorts TVB's such that the valence band maximum 
(VBM) appears as a point similarly lying along A — H . This accordingly explains why band 
touching between TVB's and BCB's occurs along A — H direction at P c . It is to be noted 
that as a direct result of H w (k), the energy gap along A — L directions is always larger than 
that along A — H directions for the whole range of P's, as shown in Fig. (2) and also the 
Supplementary Fig. S2-(e). Accordingly, neither at P c nor at any other P the band touching 
can happen along A — L directions. The effect of P is to enhance A for P < P c such that at 
the band touching points one can clearly see a rather large gap (as large as 40 meV) along 
A — L directions where H w contribution is zero. Such a warping effect and its enhancement 
at quantum critical point can be well understood by comparing Figs. 5-(b), 5(c) and 5-(d), 
in which the isocontours of energy for an arbitrary Ep fixed at 20 meV above the CBM 
(corresponding to the dashed brown lines in Fig. 2) are shown at different P's. As can be 
seen, at P c the inner and outer branches of Rashba-split conduction bands are just about 
to form the six energy pockets, whereas away from P c they form two distinct rings among 
which the outer one appears to be more significantly distorted by H w . 

Let us now address the first question, that is, why the gapless state in BiTel can be 
realized only at a certain P c but not for a range of P's. For the description of gap-closing we 
can focus on the two bands touching at P c among the four bands near the chemical potential. 
Generally, the topological phase transition in time-reversal invariant non-centrosymmetric 
systems can be described by using a two-band Hamiltonian #2x2 (k, P) = Yli=o fife, P)Tj 
where 71,2,3 are the Pauli matrices and tq is the unit matrix. The real functions /o, 1,2,3 
depend on the three momenta k = (k x , k y , k z ) and P, and are assumed to take into account 
the interaction of conduction and valence bands. The topological phase transition through 
an accidental band touching occurs if and only if the three conditions of /i,2,3(k, P) = are 
satisfied simultaneously in the (k, P) space. According to the recent work by Murakami 



35 



the band touching points, in general, form a curve in the (k, P) space because the three 
conditions /i,2,3(k, P) = cannot uniquely specify the four parameters (k, P). Therefore if 
the system is free of additional constraints other than the time-reversal symmetry, a gapless 
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metallic phase is expected to appear in a finite range of P G (P c i,P C 2) between the two 
critical pressures P c \ and P&. The occurrence of the single P c in BiTel is traced back to 
the fact that the band touching occurs along a particular direction in BZ on which the 
Hamiltonian has an additional symmetry constraint. Along the A — H line, e.g., (k x , k z ) = 
(0,7r/c), the system is invariant under the combined operation Q = TM of the T and M 
symmetries. For T = ia y K and M = ia y , Q is given by K, which imposes the following 
reality conditions: QPi2x2(k y , P)Q^ = H% x2 (k y , P) = H 2X 2(k y , P). Because of this reality 
condition, fi{k y , P) = at all points along the A — H direction. Therefore the band touching 
can be achieved when the two conditions fi^ 3 (k y ,P) = are satisfied in the (k y ,P) space. 
Since the number of conditions to be satisfied is the same as the number of parameters, a 
gapless phase appears only at a single critical point (k yjC ,P c ). 

To conclude, using the first-principles calculations we have examined the role of spin- 
orbit interaction in the giant bulk Rashba semiconductor BiTel and found that under the 
application of pressure, it leads the system to be an unusual topological insulator. The 
novel surface states and quantum critical phenomena are intriguing problems for further 
studies. The effects of the long range Coulomb interaction and disorder are, for example, 
two important issues remaining to be addressed by both the theory and experiment. 
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Methods 

Volume optimizations. To simulate the effect of P, we optimized the structure of BiTel 
crystal at various volumes ranging from V/Vq = 1 to V/Vq = 0.84, where Vq denotes the 
volume of BiTel unitcell at P am bient (o> — 4.339 A and c = 6.854 A). For a given volume, 
both the atomic positions and crystal's shape were allowed to be fully optimized until the 
magnitude of force on all ions became less than 0.005 eV/A. All the structural optimiza- 



unctional of 
. The 
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tion calculations were performed using both the LDA functional and GGA 
Perdew-Burke-Ernzerhof (GGA-PBE) as implemented in the VASP program 
corresponding Brillouin zone was sampled by a 20 x 20 x 20 /c-mesh. The respective upper 
and lower limits of bulk modulus B was estimated by fitting the LDA and GGA-PBE free 
energies E(V) to the Murnaghan equation of state 34|, using Eqn. [TJ To relate the volume 
changes to P, we then used Eqn. [2j 

Bulk electronic structure calculations. Within the same level of GGA-PBE theory, 
the electronic structures of the optimized structures were calculated using the augmented 
plane wave plus atomic orbitals (APW-LO) method as implemented in the WIEN2K pro- 



gram 



39] . For this calculations, the muffin tin radii were set to Rmt — 2.5 bohr for all 



the atoms and the maximum modulus of the reciprocal vectors K max was chosen such that 
RMTK max = 7.0. 

Surface electronic structure calculations. To calculate the surface band structure, 
for each V/Vq the corresponding APW-LO Hamiltonian was first downfolded into an effec- 
tive low en ergy 18 x 18 tight binding model using maximally localized Wannier functions 



(MLWF's) 40l442| . We chose p x , p y and p z states as the projection centers on all I, Te and 



Bi atoms to span the top 12 valence bands and the 6 lowest conduction bands by our 18 x 18 
models. Based on these realistic models, we then constructed large tight binding supercell 
Hamiltonians composed of 50 units of BiTel along its hexagonal c-axis. The surface states 
of each side is then computed by diagonalizing the whole Hamiltonian and projecting the 
corresponding eigenstates onto the MLWF's of the corresponding surface layer. 
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Fig. 1: Crystal structure, Brillouin zone and schematic diagram of band 
splitting, (a) Crystal structure of BiTel. A triple layer with Te-Bi-I is indicated by the 
purple square, (b) The relative in-plane positions of Bi, Te and I as seen along the z axis, 
(c) The corresponding Brillouin zone. The pressure-induced evolution of BiTel from a 
trivial insulator to a topological insulator is schematically drawn in (d-f). In each panel, the 
evolution of atomic p orbitals into the conduction and valence bands of BiTel is described as 
the crystal field splitting and spin-orbit coupling are turned on in sequence. As shown in (d) 
below P c , the bulk band gap Eq is trivial due to the chemical bonding. At P c , as depicted 
in (e) a gapless state is realized. Eventually, for pressures beyond P c , the band diagram in 
(f) shows a band inversion around Ep, characterizing topological insulating phase in BiTel. 
Here ±1/2 and ±3/2 indicate the ^-component of the total angular momentum. 

Fig. 2: Effect of pressure on the bulk electronic states. Electronic disper- 
sions of Rashba-split BCB's and TVB's in BiTel as hydrostatically compressed by (a) 
V/V = 1 , (b) V/V = 0.89 and (c) V/V Q = 0.86. As shown at V/V Q = 0.89 a band 
touching between BCB's and TVB's occurs along H — A direction. The dashed (brown) 
line indicates the position of a chemical potential 20 meV above the corresponding CBM 
of each case (see the related discussion). In (a) the spin directions are defined such that 
they are perpendicular to the k z well as their corresponding /c-vector. The orbital 

characteristics of BCB's and TVB's of (a), (b) and (c) are shown in (d), (e) and (f), respec- 
tively. The red, blue and green balls correspond to Bi-6p, Te-5p and l-5p states, respectively. 

Fig. 3: Effect of pressure on the surface states. Electronic band dispersions 
near Ep as obtained for the I-terminated side (top panels) and Te-terminated side (bottom 
panels) of BiTel, hydrostatically compressed by (a-b) V/Vq = 1 , (c-d) V/Vq = 0.89 and 
(e-f) V/Vq = 0.86. The Fermi surfaces corresponding to (e) and (f) are shown in (g) and 
(h), respectively. As depicted schematically in the insets, for an arbitrary Ep located at the 
middle of bulk gap, the Fermi surface at I-terminated side has a completely different shape 
from that of Te-terminated side, but interestingly for both sides, the same spin helicities 
are seen. 

Fig. 4: Volume-dependence of the free energy. Volume-dependence of free 
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energy E(V) as obtained from (a) GGA and (b) LDA calculations. To estimate the upper 
and lower limits of bulk modulus Bq and the corresponding pressure derivatives of bulk 
modulus B' Q , the corresponding GGA and LDA data are fitted to the Murnaghan equation 
of state (solid line). 

Fig. 5: Energy evolution of the conduction bands. (a) Energy isocontours 
corresponding to various chemical potentials above CBM of BiTel with V/Vq = 0.89. 
In (b), (c) and (d), the respective energy spectra of BCB's mapped into the A plane 
{k x ,k y ,k z = tt/c) of BiTel with V/V = 1, V/V = 0.89 and V/V = 0.86 are shown. The 
upper limit of energy (represented by the red color) corresponds to the Rashba energy Er of 
BCB's at the given V/Vq . The dashed lines indicate the isocontours of energy for a chemi- 
cal potential fixed at 20 meV above the CBM of each case. Arrows denote the spin directions. 
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